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Abstract : In a previous work pQ, three of us have studied a nonlinear parabolic equation 
arising in the mesoscopic modelling of concentrated suspensions of particles that are sub¬ 
jected to a given time-dependent shear rate. In the present work we extend the model to 
allow for a more physically relevant situation when the shear rate actually depends on the 
macroscopic velocity of the fluid, and as a feedback the macroscopic velocity is influenced by 
the average stress in the fluid. The geometry considered is that of a planar Couette flow. 
The mathematical system under study couples the one-dimensional heat equation and a non¬ 
linear Fokker-Planck type equation with nonhomogeneous, nonlocal and possibly degenerate, 
coefficients. We show the existence and the uniqueness of the global-in-time weak solution to 
such a system. 

1 Mechanical context and setting of the equations 

We consider here a concentrated suspension of particles in a Couette flow. Examples of such 
suspensions are numerous: tooth pastes, cements, the blood. As opposed to some other 
complex fluids such as polymeric liquids for which elaborate rheological models, based on fine 
mesoscopic physical descriptions, are available (see e.g. [SI), the modelling of concentrated 
suspensions is still in its infancy. The specific model considered here however raises interesting 
mathematical issues, mainly related to the various nonlinearities present and the coupling of 
equations at different scales. Such features are likely to be shared by a large variety of models, 
which motivates, and enlarges the scope of, the present mathematical study. 

Let us begin with some basics on the mechanical context. Depending on the concentration, 
a suspension of particles may exhibit different rheological behaviors. At low concentration, the 
suspension behaves like a newtonian fluid at rest or under weak stresses. On the other hand, 
when the suspension becomes more concentrated, the motion of each particle becomes strongly 
perturbed by the presence of the others and one observes a so-called jamming transition 
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where the sample adopts a pastelike behavior. In this transition, a macroscopic yield stress 
appears [7j. 

It is well known that when simple fluids are sheared, stress and shear rate are linked by 
a linear relation. The linear response coefficients and their relation to the microstructure of 
the fluid are well understood [H]. On the contrary, complex fluids exhibit highly nonlinear 
properties far from being understood. These nonlinear properties occur not only at high shear 
rates, where one does expect that linear response theory fails, but also at very low shear rates, 
which is more surprising. It is for instance commonly observed that for some materials (yield 
stress fluids) the shear stress a goes to a non-zero value when the shear rate goes to zero. 

In [HI, Hebraud and Lequeux proposed a model of the rheological behavior of complex 
fluids based on elementary physical processes. The system is divided in mesoscopic blocks 
whose size is large enough for the stress and strain tensors to be defined for each block. 
The size is however small compared to the characteristic length scale of the stress field. A 
mesoscopic evolution equation of the stress of each block is then written: 


i. at low shear, each particle keeps the same neighbors, and a block behaves as an Ein¬ 
stein elastic solid, in which the elasticity arises from interactions between neighboring 
particles ; 

ii. then, deformation induces local reorganization of the particles, at a given stress thresh¬ 
old Gc- Above this threshold, the block flows as an Eyring fluid : the configuration 
reached by shearing the suspension relaxes with a characteristic time Tq towards a 
completely relaxed state, where no stress is stored ; 

iii. lastly, coupling between the flow of neighboring blocks must be included. This is taken 
into account by the introduction of a diffusion term in the evolution equation, where it 
is assumed that the diffusion coefficient is proportional to the number of reorganizations 
per unit time. 


The equation proposed by Hebraud and Lequeux (HL equation in short) is written for a 
given shear rate 7 , which only depends on time: 


dtp{t,G) = -Goj{t) dap{t,G) + D{p{t)) 

](f^) . D(i 


o-j 
Dipit)) 


( 1 . 1 ) 


with 

Dipit)) = ^ [ pit, a) da. ( 1 . 2 ) 

-^0 Jjcr|>(Tc 

In the model, each block carries a given shear stress cr (cr is a real number; it is in fact an extra¬ 
diagonal term of the stress tensor in convenient coordinates). The evolution of the blocks is 
described through a probability density pit, a) which represents the distribution of stress in the 
assembly of blocks at time t. In equation (El), 1r\[-(Tc,o-c] denotes the characteristic function 
of the open set M \ [—Gc, g^ and do the Dirac delta function on M. The three terms arising 
in the right-hand side of equation ED correspond to the three physical features described 
above. When a block is submitted to the shear rate 7 , the stress of this block evolves with 
a variation rate Gq 7 where Go is an elasticity constant. When the modulus of the stress 
overcomes the critical positive value Gc, the block becomes unstable and may relax into a 
state with zero stress after a characteristic relaxation time Tq. This property is expressed by 
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the last two terms in jni). This relaxation phenomenon induces a rearrangement of the other 
blocks and this is finally modelled through the (nonlinear) diffusion term D{p{t)) d^^p. The 
diffusion coefficient D{p{t)) as given by ( 11 . 21 ) is assumed to be proportional to the density of 
blocks that relax during time Tq. The parameter a depends on the microscopic properties 
of the sample and is supposed to model the “mechanical fragility” of the material. This 
nonlinear diffusion term emphasizes the importance of collective effects in such materials. 

As mentioned above, the shear rate 7 inserted in the original HL equation depends only 
on time, and not on the space variable. It is however known from experiment that the shear 
rate in Couette flows of non-newtonian fluids is not homogeneous in space. In order to better 
describe the coupling of the macroscopic flow with the evolution of the mesostructure, we 
therefore introduce a space-dependent shear rate 7 given by the velocity gradient (which im¬ 
mediately implies that an equation of HL type holds at each point of the sample) and propose 
here the following multiscale model for planar Couette flows of concentrated suspensions (see 
Figd below) : 


pdtU{t,y) = dyT{t,y) + pd^yU{t,y) ; 

dtp{t,y,a-) = -GodyU{t,y) d^p{t,y,a) D{p{t,y)) dl„p{t,y,a) 

^ , D{p{t,y)) ^ ^ ^ 

- 7^ - P{t, y, cr) +-<50(0- ; 

To a 

Tit,y) = / ap{t,y,a)da ; 

Jr 

D{p{t,y)) = — / pit,y,a) da-, 

-^0 J|cr|>(Tc 

p > 0 ; 

U{^,y) = uoiy) , piQ,y,a) =poiy,a) 
t G (0; T), y G (0; L), u G M . 


(1.3a) 


(1.3b) 

(1.3c) 

(1.3d) 

(1.3e) 

(1.3f) 

(l-3g) 


This system is supplied with the initial condition for the probability density: 

Po > 0 , / po{y,a)da = 1 , for almost every y G ( 0 ;L) . (1-4) 

Jr 

and the no-slip boundary conditions 

U{t,0)=0, U{t,L) = V{t), for almost all t . (1-5) 

In the above equations, U{t,y) denotes the component along Cx of the velocity field (the 
flow being laminar and incompressible, the velocity field is of the form U = U{t,y) Cx), p is 
the volumic mass of the fluid and p some non-negative viscosity coefficient. The function 

V which appears in the boundary condition dni) is a continuous function on M such that 

V G L°°(M) n H[)j^(M) and H(0) = 0. The initial velocity uq lies in L^(0; L). 

The mathematical analysis of the original HL model (dl) has been the subject of pQ, where 
a more detailed presentation of the physical background and some additional references may 
be found. The main difficulties of course come from the nonlinearity in the diffusion term, 
from the presence of the singular Dirac mass as a source term, and foremost from the fact that 
the parabolic equation degenerates if the viscosity coefficient Dip) vanishes. In particular. 
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it has been shown that such a degeneracy may only occur if it is already the case at start, 
that is D[pq) = 0. And then, the situation becomes very intricate since in some particular 
cases several solutions may exist (see P). In the coupled system jni), we have to deal with 
additional difficulties due to the multiscale coupling: there are several HL-equations (roughly 
stated, one for each y) and all of them are coupled through the macroscopic equation (jl.Hal) . 
Proving the well-posedness of the Cauchy problem for the coupled system is the main purpose 
of this article. 

Before we get to the heart of the matter, we would like to comment on the diffusion term 
pdyyU in the equation of motion (II.,‘fall . Let us emphasize that this artificial viscosity has 
been added only for mathematical purposes: there is no physical reason why the fluid should 
be considered viscous. In the absence of such a regularizing term, we are unfortunately unable 
to conduct the mathematical analysis in the whole generality. 

There is however one particular situation, namely that when Uc = 0, where we are indeed 
able to study the system even without the diffusion term (i.e. with // = 0). This is the subject 
of the short Section 3. 

Notation: For given positive constants T and L, we denote H = [0; L] and = (0; T) x H. In 
the sequel, C a generic positive constant that may depend on the data but that is independent 
of t, y and a. Also to simplify the notation we shall use in the proofs the shorthands Ly, 
and La for the functional spaces LP{0;T), LP{LI), respectively. 

2 Existence and uniqueness of weak solutions 

In this section, we assume that cJc > 0, which is the physically relevant case. With an appro¬ 
priate change of scales in the coordinates and variables (see Appendix), we may equivalently 
assume that L = Tq = Uc = 1. In addition, without loss of generality, we take /U = 1 in H1.3al) . 
for, we recall, y is only here for mathematical convenience and is needed to be strictly positive. 

Let us define the velocity field 


U{t,y) = V{t)y 
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as a lifting of the boundary condition USD, and denote for T > 0 

Ut = ([0, T ]; (O)) n L 2 ([0, T ]; Hi (n)) 

and 

Vj. = C^{[ 0 ,T]-,L^{R))nL^{[ 0 ,T]-,H^{R)) . 

Setting U = u + U and denoting by V the derivative of V with respect to time, the prob¬ 
lem under consideration now reads : 


Find u G Ut and p G L“([0,T] x n L^{Q]Vt) solutions to 

/ 

p dtu - dyyU = dyT - pV{t)y\ 

< T = T{t,y)= / ap{t,y,a)da ; 

Jr 

u{0,y) = uo{y) ; 


(2.1a) 

(2.1b) 

(2.1c) 


dtp + Go {dyu + V{t)^ d„p - D{p{t,y))dl^p + ]1r\[-i,i](o-)p = 

p > 0 ; 

D{p{t,y)) = a / p{t,y,a)da] 

J\a\>l 

P(0,y,cr) =po{y,a) . 


D{p{t,y)) 


a 


6oW) ; (2.2a) 
(2.2b) 

(2.2c) 

(2.2d) 


Our main result is the following : 

Theorem 2.1 Let uq he in Lp‘{Ll) and let po satisfy the conditions 

Po > 0 ) / PoiV:'^) da = 1, for almost every y ^ Ll , 

Po £ L°°{Q xR), |(t| Po da G L^(n) , 

' Jr 


together with: 


There exists a positive constant p such that 
a inf / Po{y, a) da > p > 0 . 

yen J\a+x\>i 


xeK 


(2.3) 


(2.4) 


Then, there exists a unique global-in-time weak solution {u;p) 

u G c%R+-,L\n))nLl,(R+-Hl{n)), 

p G L°°(M+X fl;L^(M)) p e L°° {FL]Vt) , VT > 0 (2.5) 

to In addition, for such a solution, p G L(^^(M+; L°°(0 x M)) and we have 

T G L2(fl;L-(0;r))nO°([0,r];L2(Q)), 
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and 


p G C°([0,r];L2(o X M)) 


/ P{t,y,cr)da = 1, 
Jr 


for all t > 0 and y € Q, 


inf ^D{p{t,y)) > |e ^ . 
yen 

Some comments regarding the assumption (El) are immediately in order. 

Condition El obviously implies that D{pq) is bounded away from zero independently 
of y. The aim of this condition on the initial data pq is to ensure that, at any time, the 
viscosity term D{p) in H2.2al) is also bounded away from zero, so that the nonlinear parabolic 
equation satisfied by p is non-degenerate at any time (see in Lemma E21 below). 

The condition is satisfied for example when pQ is a Gaussian-like function. 

Such an assumption seems very demanding, and thus restrictive from the viewpoint of 
applications. In fact, some numerical simulations performed by one of us in |Hj show that 
even when that assumption is not satisfied at initial time t = 0 , it is indeed satisfied for 
t > 0 arbitrarily small. We are unfortunately not able to establish this fact rigorously, 
but the numerical evidence mentioned above heuristically shows that condition (El can be 
considered to be always satisfied, up to a change in the choice of the origin of times. 

The rest of this section is devoted to the proof of Theorem 12.11 The existence and unique¬ 
ness result is first proven on a small time interval with an argument based on the Banach 
fixed point Theorem. We introduce the function fFi which associates to every function u 
in ([ 0 ,T];LIq( fl)) the function r = fj^apda in L°°([0, T]; L^(n)), corresponding to the 
(unique) solution p in T] x n;L^(M)) n to H2.2I) . Then, we denote by 

fF 2 the mapping from L°°([0, T]; L^(n)) to L^([0, T]; LIq (fl)), which associates to every r in 
L°°([ 0 , r]; L^(n)), the unique solution v (in Ut ) to the heat equation : 

f pdtv- dyyV = dyT- pV{t)y on fly ; (2.6a) 

\ u( 0 ,y)=no. ( 2 . 6 b) 


We next define the mapping T on T^([0, T]; LTq (fl)) as fF = fF 2 ° ■ 

T-. L2([0,r];Fi(L!)) ^ L-([0, T]; L2(0)) ^ 

u I—> r I—> V 

and our main step consists in proving the following. 

Proposition 2.1 For every T > 0, the mapping T from L^([0, T]; Lfg (12)) into itself is well- 
defined and for T > 0 small enough, it admits a unique fixed point denoted by u. 


The proof of Proposition 12.II is organized as follows. We first check in Lemma l2.II below that 
J -2 is well-defined and that it is a Lipschitz continuous function with a Lipschitz constant that 
may be chosen arbitrarily small provided the length of the time interval is reduced (Section 
12.1 1) . In Section [2.21 we next prove that J-i is well-defined. We establish in Section [2..11 that 
J-i is a Lipschitz continuous function with a locally bounded Lipschitz constant with respect 
to time interval. Therefore the composed mapping T is contracting on small enough time 
interval. The existence and uniqueness of a solution on a small time interval follows by the 
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Banach fixed point theorem. Finally in Section 12.41 we deduce the existence and uniqueness 
of the global-in-time solution. 

Henceforth, and unless otherwise stated, the initial condition pQ is hxed and it satisfies 
the assumptions of the statement of the Theorem. 

2.1 The map T2 is a contraction on [0,T] for T small enongh 

Lemma 2.1 For every T > 0, the mapping T2 is Lipschitz continuous from L°° ([0, T]; 
to L^([0, T]; Hg (H)) , and the Lipschitz constant goes to 0 with T. 

Proof of Lemma l2.lt We hrst observe that the mapping T2 is well-dehned. Indeed, for every 
function r in L°°([0, T]; L^(H)), dyT G L°°{[ 0 ,T]-, and therefore, the existence and 

the uniqueness of a solution v G Ut of the heat equation (EH) is a standard result. Let now ri 
and T 2 be two functions in L°°([0, T]; L^(H)), and let us denote vi = F 2 {ti) and V 2 = F 2 {t 2 )- 
We also set v = vi — V2 and r = ri — r2 . Then, v satisfies 

pdtv— dyyV = dyT on LIt ; (2.8a) 

< u(0 ,y)=0 ; (2.8b) 

v{t,0) = v{t,l) = 0 , (2.8c) 

and if we apply Equation (EH to V and integrate over H we get 

+ = (2.9) 

By the Cauchy-Schwarz and the Young inequalities, we obtain for t G [0;r], 

P f + f ( f \dyv\‘^ dy) ds < [ ( [ dy) ds , 

Jn Jo ^Jn ' Jo ^Jn ' 

and therefore by the Poincare inequality 

ll'J^llL2([o,T];//i(n)) < 2\/r||T||ioo([o^'r].i2(Q)) . (2.10) 

0 


2.2 The map T\ is well defined 

Equation (EH) with the variable y frozen has been studied in pQ. For the sake of consistency 
we now recall : 

Proposition 2.2 |7[ Theorem 1.1] (Global-in-time existence for all y) 

For almost every y in Q, let b{-,y) be a given function in L^q^(M+), and let po such that : 

Po{y,-) e , po(y,-)>0, fpo{y,a)da = l and f \a\po da <+00 , 

( 2 . 11 ) 
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and 


D{po{y)) > 0 . (2.12) 

Then, for every T > 0 and for almost every y in 11, there exists a unique solution p = p{t, y, a) 
in L°°([0, T]; L^(M) n L°°(M)) n L^([0, T]; (M)) to the equation 

dtP =-b{t,y) d„p + D{p{t,y))dl„p- lR\[-i,i]icr)p + 6o{cr) ; (2.13a) 

p > 0] (2.13b) 

D{p{t,y)) = a [ p{t,y,a)da; (2.13c) 

J |o-|>l 

p(0,y,cT) =po(y,o-) • (2.13d) 


In addition, for almost every y in Q, we have 

• / P{t^ Vi O') da = 1 for all t > 0 , 

Jr 

• for all T > 0, 


max ||p(t,y,-)||L“ < ||j>o(2/, OIIl- + 

\/vr 


. p(-,y)eCO([0,r];Li(M)nL2(M)), 

. D{pi;y))€C^{[0,T]), 

• for every T > 0 there exists a positive constant r]{T, y) sueh that 

^^I^^D{p{t,y)) > v{T,y) . 


(2.14) 


(2.15) 


• for almost all y, {t,a) i— > ap{t,y,a) G L°°([0,T];L^(M)), so that the average stress 
T{-,y) is well defined by U.dd) m L-(M+). 

We now fix some initial condition po satisfying the conditions (tOl and (|2.4I) (thus a fortiori 
the conditions (tTTTl) and (tTT^ i and set 

b{t, y) = Go {dyu{t, y) + V(t)) 

for u G L^([0, T]; (H)). In view of Proposition 1221 we know the existence and uniqueness 

of a solution p to (|1.3b|) for given u. Our next step now consists in analyzing the dependence 
on the variable y. 

Lemma 2.2 (Uniform-in-y a priori estimates on p) Let T > 0 be given. We assume 
that the initial data po satisfies (tO) and 

inf D{po{y)) > 0 . 
y&n 

Notice that \2.4\) is not needed, but flTil) and h2.1^) ) are fulfilled. Then, if we denote by p 
the unique solution to (TOl) given by Proposition m we have : 
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(i) p E L°°([0,r] X n^LlnL^) with 


l|p||i,°°([o,r]xnxR) < l|Po||L°°(nxR) + 


Vt 


TT 


and 


p{t, y, a)da = 1 for all t > 0, for almost every y in Q. 


(2.16) 

(2.17) 


(a) The stress T is in L'^{Q.,L°°{[Q,T])) (hence in L°°{[{),T]]L'^{Q))). 

(Hi) If in addition po satisfies the non-degeneracy condition (EH), we have 

inf D{p{t,y)) > - e~'^ n , 
o<t<T - 2 ' ’ 

yen 

and 


(2.18) 


sup^ (||po||L-(nxR)(^+ 7") +• (2-19) 

Proof of Lemma 12.21 

To prove Assertion (i), we use the estimates obtained in pQ with the variable y kept frozen. 
The assumptions on pQ ensure that p is in for almost every y, that (ItTTI) holds, and that 
(ITTH) holds by virtue of E Proposition 1.1, Eq.(1.8)]. Estimate (j2.16j) follows. 

Assertion (ii) follows from ^ Proposition 1.1, Eq.(1.9)] : for almost every y in fl, 

sup f \a\pda < f \a\poda-I VT + \\b\\^A , 


with b = b{t, y) = Go {dyu(t, y) + V(t)) . Then 


sup / \a\pda < \\apo\\ l^lI) + I|P||l“([o,t])) + VTGo\\dyu\\L 2 (^^^) , 

o<t<T Jm. ^ 


with 


G(r,a,||E|Ucc([o,T])) = G'or 


v/T 

L-([0,T]) + - 7 ^ 



Eor Assertion (hi), following E Proof of Lemma 3.1], we define 


t*(y) = inf i t > 0; / po{y,a)da = 0 

./ cr+frfb(s,i/)ds >1 


ct+/q b(s,j/)ds|>l 

Because po satisfies the non-degeneracy condition (EH), we have for all y, t*{y) > 0 and the 
support of Po is contained in the interval —1 — b(s, y) ds, 1 — b{s, y) ds . 
Moreover, for any T = T{y) < 


min ^(^(^, 1 /)) > min / po{y,a)da. (2.20) 

0 <t<T 2 ^—^—^J\a+f(b(s,y)ds\>l 
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The assumption (tmi on pq ensures that t*{y) = +oo for almost every y in O. Therefore 
holds true on any time interval T > 0 independently of y and (ITTKl) is an immediate 
consequence of by using (|23I). 

Finally (ITT^ follows in a very standard way from (ITTKl) and ^ Equation (3.7)], multi¬ 
plying (ll.3bll by p, integrating over M with respect to a, and using H2.18I) and the previous 
bounds. <0 


2.3 The map J^i is Lipschitz continuous 

Lemma 2.3 For every T > 0, the mapping is Lipschitz continuous from L'^ (^[0,T]‘, HQ{n)^ 
to L°°([0,T];L^(n)), and the Lipschitz constant is a locally bounded function of T. 


Proof of Lemma lot Let ui and U 2 be two functions in L'^ [[0,T]-, , and let ri = 

Fi{ui) and T 2 = Fi{u 2 )- We denote by pi, i = 1, 2, the unique solution to (12.21) corresponding 
to Ui whose existence is guaranteed by ProDosition l2.2l and Lemma lT^ We also set v = Ui—U 2 , 
q = Pi—P 2 and r = ri — T 2 . Recall that, for i = 1, 2, r* = api da. We formally multiply 
equation (iO) by a and integrate it over M with respect to cr to find 


dtn + Ti = Go {dyUi + V{t)) + a Pi da 

J |(t |<1 

T(0,y) = / apoda . 


The argument may be made rigorous with the help of a standard cut-off argument as in 
Subtracting the equations satisfied by ri and T 2 yields 


dtT + T = Go dyV + aqda ] 

r(0,y) = 0 . 

We then apply r to and integrate over to obtain 

f + f = ^0 f + f ( f crqda] dy . 

^ntjQ jq jq jq y./|cr|<i y 

Using the Young and the Cauchy-Schwarz inequalities, we have 


1 


Go, 


\9yvr\ < + —\\dyv\\^2^^) 


and 


/ r / aqda \ dy 
In / 


thus 


fL 

dt 


^ / \Q\‘^da) ' jrj dy 

[ |rp < G^\\dyv\\l2^n) + / (/ 

J J J 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 
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and 


(2.24) 


oSup^ll'^lli^cn) < GiWdyvWl^^^^) + hWh^arxi 


Let us now admit for a while that 

ll^lli2(nyxR) — G{T)\\dyv\\'j^2(^Q^^, (2.25) 

with C{T) being a locally bounded function of T. Inserting (12.251) into (|2.24l) . we obtain 

sup ||'r||i2(Q) < {Gq + C'(T))||cIj,u||^2(q ) , (2.26) 

and therefore, the mapping J^i is indeed Lipschitz continuous. 

In order to establish (12.2511 . we subtract the equations (1^ satisfied by pi and p 2 respec¬ 
tively to deduce that q solves 


dtq = - Go dyUi d^q - Gq dyV d„p 2 - GqV (t) d^q + D{q) dl„pi 

+D{p 2 ) dl^q - q + 5o{a) ; 

9 ( 0 , y,o-) = 0 , 


(2.27a) 

(2.27b) 


for almost every y in 12. Then, we apply (ITTTl) to q, and integrate with respect to a, to obtain 

+ o{p2) f 

^ Ot J-R J\a\>l JR 

= GodyV [ P2daq-D{q) f d^pi d^q + ^^^q{t,y,0) . (2.28) 

Jr Jr a 

By the Cauchy-Schwarz and the Young inequalities and using the bound from below (ItTsI) 
on D{p 2 ), we have 


dyV / P2d„q 

Jr 


< \dyv\ \\P2 \\lI Wdaqhl 

\\P2\\l2 \dyv(^ + 


< 


3^0 „2 ,2 , D{p2) 


< 


4D{p2) 
3 Go 

2 rj 


3 Go 
y/aT 


\\d.q\\ 


2 

LI 


(IIpoIIl^^ H-, (2.29) 


thanks to the L°° bound (ITTHl) on p 2 and the fact that / P 2 {t, y,a) da = 1. In a similar way, 

Jr 

\D{q) [ daPidaq\ < \D{q)\ WdaPihl Wdaqhl 

Jr 


we obtain 


3 Tit-),' m 2 no \\2 , G>{p 2 ) ,,2 


< e^ \Diq)\'^ \\daPi\\l2 + 

Ip 6 




LI ■ 


(2.30) 


As / q = 0 (recall / pi = / P 2 = 1), we may write 

t/M. t/]R. «/M 


|L>(g)| = a 


/ ’ 

= a 

/ ’ 

J crj>l 




(2.31) 
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Thus, inserting (lOTl) into (12201), we obtain 


\D{q) [ daPid^ql < ^ || ^ ||^2 || 9<tPi 11^2 + \\daq\\l 2 . (2.32) 

Jm. V 3 

On the other hand, from the Sobolev embedding L°°(M), we know 

\\<lh~<^{\\q\\l.+\Mhji , (2.33) 

and next using successively Young’s inequality, (lOTl) . (ITTHl) again and the fact that D{p 2 ) < 
a, we find 


i-\D{q)\\q{t,y,0)\ < WqhliWQWh Hdaqllh)^ 


a 


< 


( D{P2) 

V 3 


+ 


4T»(p2) 


) Ik lli^ + 


D{P2) 


\\daq 


|2 

\lI 


(a 


- \3^ 2r]^ ) 

Inserting (121201), (t^^ and (|nni) in (12I2S1), we have 


DiP2) 


Wdqrql 


Ll 


(2.34) 


I 


2 ^Gn T 


VaT 




\dyV\ 


2r] 

( Ol 3 O UC T'iir\ ii9 

w + 73- e +- d^pi I2 

6 Irj rj 


3 


q Wli 


for all t in [0; T] and almost every y in 0. Applying the Gronwall lemma, we obtain 

3Gr? 


/ q‘^{t,y,cr) 
Jr 


X exp 


da < (ikolU-, + 

2aT 3„t’ 6 a^ 


J \dyv\'^ dt 


+ -Te^ + 

p 


j II daPi 11^2 dt^ , (2.35) 


for almost every {t,y) in O-r. We now integrate the above inequality over and we use the 
bound (12.191) on d^pi in L^(L|,^) to deduce (|2.25l) with 


C{T) = 


T 


V 

X exp 


Te‘ (||p„||,., + 2l^) X 


2aT 3 


6 


—+ -re^ +- e^C 

6 p p 


(2.36) 


with an explicit expression for C'{T) coming from the right-hand side of (12.191) . namely : 

c'(T) = Y-(||p„|U„(l + r) + ^r3/^) , 


0 
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2.4 Global-in-time existence 


Let us assume that there exists some finite t* such that the system admits a solution {u*\p*) 
on [0,t*[ that ceases to exist after the time t* (at least in the appropriate functional spaces 
prescribed by our notion of solution). 

Let 0 < to < We consider the following Cauchy problem starting at time to : 


p dtu - dyyU = dyT - pV{t) y ; 

(2.37a) 

T(t,y) = / o-pdcr ; 

(2.37b) 

J M. 

n( 0 , y) = uo{y) := u*{to, y) on 0 ; 

(2.37c) 

u{t, 0 ) = 0 , u{t, 1 ) = 0 , 

(2.37d) 


which is coupled to 

dtp + Go {dyu + V{t)) d„p - D{p{t,y)) dl„p + ]1k\[-i,i](u)p = (5o(cr) ; (2.38a) 

< p > 0 ; (2.38b) 

P(0, y, cr) = po{y, a) := p*{to, y, a) . (2.38c) 

As a first step we prove that the new initial condition {uo;po) satisfies the assumptions in 
Theorem 12.1 1 The only point to be checked is that po fulfills (1221). Indeed, uo is in L^(n) 
and po satisfies dH, thanks to the bounds in Lemma 1221 which hold true for p* on [0,t*). 
Actually we prove the more general : 

Lemma 2.4 Let (tt*;p*) be the solution to the coupled system on [0; t*) under the assumptions 
of Theorem, \2. /I for the initial data. Then, for every 0 < T < t*, we have 

inf f p*{t,y,a)da > ^e~'^ . (2.39) 

y|cr+Y|>l 2 O 

XGR 


Proof of Lemma 12.4t The proof follows from a comparison principle and is inspired 
from pp. It is reproduced here for the reader’s convenience. We denote by p_ the solution to 
the linear equation : 


f dtp- = 

1 P-{0,y,cr) = 

It is well-known that p_ 

P-{t,y,cr) = e~^ 

with 


- Go {dyU + V (t)) daP- + D{p*{t, y)) dl„p- - p_ ; 
Po{y,cr) ■ 

is given by 


r+oo 


Po(y, a') <p / .. . (o- -cf' - ({t, y)) da' 

Jo D(j}*{s,y)}ds 


Tu{x) = 


1 


V 

To = (5o , 


exp ( - 7 ^) if z/ > 0 ; 


(2.40) 


(2.41) 
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and 


^{t,y) = Go [ {dyu{s,y) + V{s)) ds . 

Jo 


Since p- < p*, by the maximum principle, we get the bound from below 


'k-xl>i 

> f 


p*{t,y,a)da 

p-{t,y,a)da 

J\a-x\>l 

> e~^ / Po{y,o'')( / p /—-t - (a — a'— ^(t,y)) da) da , (2.42) 


for every X in M. As in pQ , we introduce the interval = [ —1 — + x, 1 — + x]- 


The function a ^ ip 


-{a — a' — i{t,y)) is a Gaussian probability density with 


\/2 fo Dip*i^>y))ds 

mean a' + ^{t, y) and squared width 2 f* D(p*(s, y)) ds. Therefore, for every u' e M \ 
we have 

/ ‘P /—TT -(<7 — a' — f(t, v)) da > - , 

-^0 2 

which implies 


(|2.42D >l-e ^ f po{y,a')da' = Je f 
^ Jr\Kc., ^ j\a' 


<R\K^,^ 

And we conclude using (ini) . 




Po{y,a')da' . 


0 


Completion of the Proof of Theorem 12.11 

In view of Lemma im we may apply the Banach fixed point theorem as in the proof of 
Proposition 12 .II and deduce the existence of a unique solution to the Cauchy problem 
(lOHl) on the time interval [toi^o + for some small enough k > 0. Of course, this solution 
coincides with (u*;p*) by uniqueness. We now show that k may be chosen independently of to 
in (0;t*). Therefore the solution exists beyond the time t*, which contradicts the finiteness of 
t*. The constant k depends on the Lipschitz constant for the mapping T. Because of (10)1) 
the Lipschitz constant of the mapping T 2 is clearly independent of the initial time to and can 
be made arbitrarily small using (10)1) provided the length of the time interval is taken small. 
We thus focus on the Lipschitz constant of J^i, and now show it is bounded uniformly in t*. 
Thus, the condition on k such that .F = .T 2 o .Fi is a contraction on [to; to + k] is independent 
of to in (0;t*), which concludes the proof. 

We revisit carefully the proof of Lemma 12.41 which is the crucial step for checking the 
assumptions of the Banach fixed point theorem on small time interval. We go back to the 
proof of (12.2511 . The only modifications are in the proofs of estimates dHH), (IT^T^ and 
(|2.44|) as follows. In view of the uniform estimate given by Lemma 1231 the quantity D{p*) is 
bounded from below by y exp(—t*) /2 and the L°° and norms of p* in the sense of (12.151) 
and (1tT91) are bounded uniformly in terms of t*. (In all these bounds T may obviously be 
bounded by t*.) Therefore the Lipschitz constant of J-i given by (12:2Hl) is bounded uniformly 
in t*. 0 
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3 The case (Jc = 0 


In the situation examined so far, that is when Uc > 0, we have only been able to show the 
well-posedness of the coupled system when the macroscopic equation has a positive diffusion 
coefficient /x. This is a mathematical artefact, apparently related to our technique of proof. 
Our aim in the present section, as announced in the introduction, is to mention that the 
coupled system is also well-posed in the particular case when /x > 0 and cic = 0 . 

We again scale out the variables y and f, together with the function u as explained in 
Appendix, which amounts to taking Tq = L = 1. In the present case when Uc = 0, we 
note that, for a given b{t,y) £ L^(0)), the unique solution of (I2.13|) provided by 

Proposition reads 

Pit,y,cr) = e~^ [ pQ{a')Qat{cr - o' -x{t,y)) da' (3.1) 

+ / Qa{t-s) (<7 - x(T y) + x{s, y)) ds, 

Jo 

where x(t, y) = / b{s,y)ds and where {Qt)t>o denotes the heat kernel 

Jo 

Qo{o) = dQ{a) . 

If we multiply eu by a and integrate over the real line, we obtain 

dtT{t,y)+T{t,y) = b{t,y) . (3.2) 

For Oc = 0, the multiscale Hebraud-Lequeux model is therefore equivalent to the so-called 
Maxwell model [HI, and the coupled system under consideration reads 

pdtu = pL dyyU + dyT - pV (t) y ; 

< dtT + T = Gq dyU + GqV( t) ; (3.3) 


^ n(-,0) = xx(-, 1) = 0 . 

The latter is a linear system for which it is easy to prove global existence and uniqueness 
in convenient functional spaces such that dyU £ L^q^(M'’'; L^(n)), whatever y > 0. Global 
existence and uniqueness for the multiscale Hebraud-Lequeux model immediately follows. 


Appendix: Non-dimensionalized equations 


We first scale the space and time variables in order to work with dimensionless constants and 
with a reduced number of parameters. We introduce the new dimensionless variables 


t' 



y 

L’ 


/ 


a 


a 

? 

Oc 
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and the dimensionless rescaled functions 


T ^0 j j / / ^ f f f ^ f 

U = — U , p = cTcP , T = — = j a p da 


together with the corresponding dimensionless parameters 


J _ P ^ ^ _ ^0 

P _T-i95 ^ 2 ’ ^0 


/ p 
p = 


Z i 9 5 9 5 ^ u : rri 

0 cri CTc Toac 

Note that p' is actually the so-called Reynolds number. We also define 

D'{p') = a' f p'da' . 

J\a'\>l 

Then, equations (Oa and (lObll respectively read : 

pdtU' - pdyyll' = dyr' 


(3.4) 


and 

dtp' = -Go dyU' d„p' + D{p') dl„p' - ]lR\[-i,i](o-)y -b 6o{a) . (3.5) 

Of course the corresponding change of scales and variables are also applied to the initial 
conditions uq and po. The new function V entering in the boundary conditions of \J' has to 
be changed according to 

All the primes are omitted in the body of the article in order to lighten the notation. 
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